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This application claims the benefit of United States Provisional Patent Application No. 
60/524,326, filed on November 20, 2003, the disclosure of which is hereby incorporated by 
reference. 

Portions of the disclosure of this patent document contain material that is subject to 
5 copyright protection. The copyright owner has no objection to the facsimile reproduction 
by anyone of the patent document or the patent disclosure as it appears in the Patent and 
Trademark Office file or records, but otherwise reserves all copyright rights whatsoever. 

BACKGROUND OF THE INVENTION 

1 . FIELD OF THE INVENTION 

1 0 The present invention relates to the implementation of a method and apparatus for 
estimating battery charge power and discharge power. 

2. BACKGROUND ART 

A number of high-performance battery applications require precise real-time estimates 
of the power available to be sourced by the battery pack. For example, in Hybrid Elec- 

1 5 trie Vehicles (HEVs) and Battery Electric Vehicles (BEVs), the vehicle controller requires 
continuous up-to-date information from the Battery Management System (BMS) regarding 
the power that may be supplied to the electric motor from the battery pack, and power 
that may be supplied to the pack via regenerative braking or by active recharging via the 
motor. One current technique in the art, called the HPPC (Hybrid Pulse Power Character- 

20 ization) method, performs this task of estimation by using the voltage limits to calculate 
the maximum charge and discharge limits. As described in the PNGV (Partnership for 
New Generation Vehicles) Battery Test Manual, Revision 3, February 2001, published by 
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the Idaho National Engineering and Environmental Laboratory of the U.S. Department of 
Energy, the HPPC method estimates maximum cell power by considering only operational 
design limits on voltage. It does not consider design limits on current, power, or the battery 
state-of-charge (SOC). Also the method produces a crude prediction for horizon At. Each 
5 cell in the battery pack is modeled by the approximate relationship 



v k (t) = OCV{z k (t))-Rxi k {t), (1) 

where OCV (z k (t)) is the open-circuit- voltage of cell k at its present state-of-charge z k (t) 
and R is a constant representing the cell's internal resistance. Different values of R may 
be used for charge and discharge currents, if desired, and are denoted as R^ z and R dis , 
1 0 respectively. 

Since the design limits v m \ n < v k (t) < v mBX must be enforced, the maximum discharge 
current may be calculated as constrained by voltage, as shown below 

•dis,volt 

ocv (*(«)) 

The maximum magnitude charge current may be similarly calculated based on voltage. 
1 5 Note, however, that charge current is assumed negative in sign by convention employed in the 
present invention (although the opposite convention may be used with minor modifications 
to the method), so that maximum-magnitude current is a minimum in the signed sense. It is 

•chg,volt _ 

OCV (**(*)) 
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Pack power is then calculated as 



This prior art charge calculation method is limited in several respects. First, as noted 
5 above, the method does not use operational design limits on SOC, maximum current, or 
maximum power in the computation. More importantly, the cell model used is too primitive 
to give precise results. Overly optimistic or pessimistic values could be generated, either 
posing a safety or battery-health hazard or causing inefficient battery use. 

What is desired is a new method and appartus for battery charge estimation based on a 
1 0 better cell model. Such a cell model would be combined with a maximum-power algorithm 
that uses the cell model to give better power prediction. The new method would also take 
in operational design limits such as SOC, current, and power. 
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SUMMARY OF THE INVENTION 

The present invention relates to a method and an apparatus for estimating discharge 
and charge power of battery applications, including battery packs used in Hybrid Electric 
Vehicles (HEV) and Electric Vehicles (EV). One embodiment is a charge prediction method 
5 that incorporates voltage, state-of-charge, power, and current design constraints, works for a 
user-specified prediction horizon At, and is more robust and precise than the state of the art. 
The embodiment has the option of allowing different modeling parameters during battery 
operation to accommodate highly dynamic batteries used in Hybrid Electric Vehicles (HEV) 
and Electric Vehicles (EV) where such previous implementations were difficult. 
1 0 An embodiment of the present invention calculates maximum charge/discharge power by 
calculating the maximum charge/discharge current using any combination of four primary 
limits: 

1. state-of-charge (SOC) limits 

2. voltage limits 
1 5 3. current limits 

4. power limits 

In one embodiment, the minimum absolute charge/discharge current value from the cal- 
culations using state-of-charge (SOC), voltage, and current limits is then chosen to obtain 
the maximum absolute charge/discharge power. In one embodiment, the maximum absolute 

20 charge/discharge power is checked to ensure it is within the power limits. In one embodi- 
ment, the maximum absolute charge/discharge power is calculated in a way as to not violate 
any combination of the limits that may be used. 

Prior methods do not use SOC limits in their estimation of maximum charge/discharge 
power. The present invention incorporates the SOC of the battery cell or battery pack to 

25 estimate the maximum charge/discharge current. The estimation explicitly includes a user- 
defined time horizon At. In one embodiment, the SOC is obtained by using a Kalman filter. 
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The SOC that is produced by Kalman filtering also yields an estimate of the uncertainty 
value, which can be used in the maximum charge/discharge calculation to yield a confidence 
level of the maximum charge/discharge current estimate. 

Methods of the present invention improve upon prior art estimation of power based on 
5 voltage limits. In the present invention, voltage limits are used to calculate the maximum 
charge/discharge current in a way that includes a user-defined time horizon At. Two pri- 
mary cell model embodiments are in the present invention for the calculation of maximum 
charge/discharge power based on voltage limits. The first is a simple cell model that uses a 
Taylor-series expansion to linearize the equation involved. The second is a more complex 

1 0 and accurate cell model that models cell dynamics in discrete-time state-space form. The 
cell model can incorporate a variety of inputs such as temperature, resistance, capacity, etc. 
One advantage of using model-based approach is that the same cell model may be used in 
both Kalman filtering to produce the SOC and the estimation of maximum charge/discharge 
current based on voltage limits. 

1 5 Embodiments of the present invention also include methods of charge estimation based 
on any combination of the voltage, current, power, or SOC limits described above. For 
example, charge estimation can be based on voltage limits only, or combined with current 
limits, SOC limits and/or power limits. 

Embodiments of the present invention are directed to a power estimating apparatus that 

20 takes in data measurements from the battery such as current, voltage, temperature, and 
feeding such measurements to an arithmetic circuit, which includes calculation means that 
performs the calculation methods disclosed in the present invention to estimate the absolute 
maximum charge or discharge power. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

These and other features, aspects and advantages of the present invention will become 
better understood with regard to the following description, appended claims and accompa- 
nying drawings where: 

5 Fig. 1 A is a flow chart that outlines the maximum discharge estimation according to an 
embodiment of the present invention; 

Fig. IB is a flow chart that outlines the minimum charge estimation according to an 
embodiment of the present invention; 

Fig. 2 is a schematic block diagram showing the sensor components of a power esti- 
1 0 mating embodiment of the present invention; 

Fig. 3 is an example plot of open-circuit- voltage (OCV) as a function of state-of-charge 
for a particular cell electrochemistry; 

Fig. 4 is an example plot showing the derivative of OCV as a function of state-of-charge 
for a particular cell electrochemistry; 
1 5 Fig. 5 is a plot showing the voltage prediction using the cell model of the present 
invention; 

Fig. 6 is a zoom-in of the plot of voltage prediction for one UDDS cycle at around 50% 
state-of-charge; 

Fig. 7 is a state-of-charge trace for cell test; 
20 Fig. 8 is a plot comparing static maximum power calculations as functions of SOC for 
the PNGV HPPC method and Method I of the present invention; 

Fig. 9 is a plot showing that discharge power capability estimates for cell cycle test 
comprising sixteen UDDS cycles over an SOC range of 90% down to 10%; 

Fig. 10 is zoomed-in plot of Fig. 9, showing about one UDDS cycle; 
25 Fig. 1 1 is a plot showing charging power capability estimates for cell cycle test com- 
prising sixteen UDDS cycles over an SOC range of 90% down to 10%; and 
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Fig. 12 is zoomed-in plot of Fig. 11, showing about one UDDS cycle. 



DETAILED DESCRIPTION OF THE INVENTION 



Embodiments of the present invention relate to battery charge estimation for any battery- 
powered application. In one embodiment, the estimator method and apparatus find the 
maximum absolute battery charge and/or discharge power (based on present battery pack 
5 conditions) that may be maintained for At seconds without violating pre-set limits on cell 
voltage, state-of-charge, power, or current. 

Figs. 1A and IB illustrate an overview of the embodiments of the present invention. 
Fig. 1 A shows a method for finding the maximum discharge power for a user-defined time 
horizon At, i.e. how much power may be drawn from the battery continuously for use for 

1 0 the next At time period. In vehicle applications, accurate estimation of maximum discharge 
power can help prevent the hazardous occurrence of over-drawing the battery. 

In step 10, the maximum discharge current is calculated based on pre-set limits on 
state-of-charge. The estimation explicitly includes a user-defined time horizon At In one 
embodiment, the SOC is obtained by using a Kalman filtering method. The SOC that is 

1 5 produced by Kalman filtering also yields an estimate of the uncertainty value, which can 
be used in the maximum charge/discharge calculation to yield a confidence level of the 
maximum charge/discharge current estimation. In another embodiment, a simple state-of- 
charge is used. Step 10 is further described in the section titled "Calculation Based on 
State-of-Charge (SOC) Limits;' 

20 The maximum discharge current is calculated based on pre-set limits on voltage in 
step 12. The present invention has two main cell model embodiments for the calculation 
of maximum charge/discharge power based on voltage limits, although it is understood 
that other cell models could be used. Both overcome the limitation of prior art discharge 
estimation methods of giving a crude prediction of time horizon At . The first is a simple cell 

25 model that uses a Taylor-series expansion to linearize the equation involved. The second is a 
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more complex and accurate cell model that models cell dynamics in discrete-time state-space 
form. The cell model can incorporate a variety of inputs such as temperature, resistance, 
capacity, etc. The two cell models are further described in the section titled "Calculation 
Based on Voltage Limits." 
5 Then in step 14, the maximum discharge current is calculated based on pre-set limits on 
current. In step 16, the minimum of the three calculated current values from steps 10, 12, and 
14 is chosen. It is understood that the execution order of steps 10, 12, 14 is interchangeable. 
It is further understood that any combination of steps 10, 12, and 14 may be omitted, if 
desired, in an implementation. Using the chosen discharge current value, step 18 calculates 

1 0 the maximum discharge power. The calculated pack power may be further refined in order 
to not violate individual cell or battery pack power design limits. 

Fig. IB shows a method for finding the maximum absolute charge power for a user- 
defined time horizon At , i.e. how much power can be put back into the battery continuously 
for the next At time period. The details and progression of the method mirror that of 

1 5 Fig. 1 A. Since charge current is considered to have a negative sign, the maximum absolute 
current is the minimum current in a signed sense. In step 20, the minimum charge current 
is calculated based on pre-set limits on state-of-charge. Again the SOC can be a simple one 
or one obtained using the Kalman filtering method. Then the minimum charge current is 
calculated based on pre-set limits on voltage in step 22 in accordance with a cell model, 

20 such as one of the two cell models described in the present disclosure. Then in step 24, the 
minimum charge current is calculated based on pre-set limits on current. Then, in step 26, 
the maximum of the three calculated current values from steps 20, 22, and 24 is chosen. 
Note again that the execution order of steps 20, 22, 24 is interchangeable. It is further 
understood that any combination of steps 20, 22, 24 may be used and any of the steps 20, 

25 22, 24 may be omitted, if desires, in an implementation. Using the chosen charge current 
value, step 28 calculates the minimum charge power. The calculated pack power may be 
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further refined in order to not violate individual cell or battery pack power design limits. 

It is noted that modifications may be made to the method embodiments as shown in Figs. 
1 A and IB. For example, any or all of the current calculation steps based on state-of-charge 
and voltage limits may be removed. Also, the present invention discloses several methods 
5 of calculating maximum absolute charge and discharge current based on state-of-charge, 
voltage limits, and current limits. 

One embodiment of the present invention estimates the maximum absolute charge and/or 
discharge power of a battery pack. The battery pack may be, for example, a battery pack 
used in a hybrid electric vehicle or an electric vehicle. The embodiment makes a number 
1 0 of denotations and limits, including: 

- using n to denote the number of cells in the target battery pack where an estimation of 
charge and/or discharge power is desired; 

- using v k (t) to denote the cell voltage for cell number k in the pack, which has operational 
design limits so that v m \ n < Vk(t) < u max m ^st be enforced for all k : 1 < k < n; 

15 - using z k (t) to denote the state-of-charge for cell number k in the pack, which has 
operational design limits z m \ n < z k (i) < z max that must be enforced for all k : 1 < k < n\ 

- using pk(t) to denote the cell power, which has operational design limits so that 
Pmin < Pk(t) < Pmax must be enforced for all k : 1 < k < n\ and 

- using i k (t) to denote the cell current, which has operational design limits so that 
20 i min < i k (t) < imax must be enforced for all k : 1 < k < n. 

Modifications can be made in alternate embodiments. For example, any particular limit 
may be removed if desired by replacing its value by ±oo, as appropriate. As an another 
example, limits such as v mdXi v m[Q , * max , z m[Qy w, i m \ n , Pmax, Pmin may furthermore be 
functions of temperature and other factors pertaining to the present battery pack operating 
25 condition. In one embodiment, it is assumed that the discharge current and power have 
positive sign and the charge current and power have negative sign. Those skilled in the 
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art will recognize that other sign conventions may be used, and that the description of the 
present invention can be adapted to these conventions in a forthright manner. 

In one embodiment, the model used for predicting charge assumes that the battery 
pack comprises n 8 cell modules connected in series, where each cell module comprises n p 
5 individual cells connected in parallel and n 3 > 1, n p > 1. Other configurations are possible 
and are accommodated by slight modifications to the method as described. 

FIG. 2 is a schematic block diagram showing the sensor components of an embodiment 
of the present invention. Battery 40 is connected to load circuit 48. For example, load circuit 
48 could be a motor in an Electric Vehicle (E V) or Hybrid Electric Vehicle (HE V). In some 

1 0 embodiments circuit 48 is a circuit that provides power and/or draws power. Measurements 
of battery and individual cell voltages are made with voltmeter(s) 44. Measurements of 
battery current are made with ammeter 42. Battery and individual cell temperatures are 
measured by temperature sensor(s) 46. Voltage, current and temperature measurements 
are processed with arithmetic circuit 50. Arithmetic circuit (estimator means) 50, takes in 

1 5 the measurements from the sensor components and perform the calculation methods of the 
present invention for power estimation. In some embodiments, temperature is not needed 
in the calculation methods. 

1 Calculation Based on State-of-Charge (SOC) Limits 

As shown in steps 10 and 20 of Figs. 1A and IB, embodiments of the present invention 
20 calculate the maximum charge/discharge current values using SOC limits. Various embod- 
iments also have the explicit inclusion of a time horizon At in the calculation. The SOC 
limits are included as follows. First, for a constant current i k , the SOC recurrent relationship 
is described as: 

z k (t + At) = z k (t) - (viAt/C)i kj (4) 
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where z k (i) is the present SOC for cell k, z k {t + At) is the predicted SOC At seconds 
into the future, C is the cell capacity in ampere-seconds, and rji is the Coulombic efficiency 
factor at current level i k . Here, for simplicity of presentation, it is assumed that m = 1 for 
discharge currents and = 77 < 1 for charge currents. 
5 If there are design limits on SOC such that z m \ n < z k (t) < z max for all cells in the pack, 
then current i k can be computed such that these limits are not exceeded. Simple algebra 
gives limits based on the SOC of each cell: 

■dis,soc _ z k(t) ~ z m\u 



At/C 

•chg,soc 

r)At/C ' 



Vin> = „a,/^ • (6) 



1 0 The pack maximum absolute currents — based only on cell SOC — are then 



C = "in (£5) 

•chg,soc _ / .chg,soc\ 



This method assumes that there is a valid SOC estimate available for every cell in the pack. 
If this is not the case, then an approximate remedy would be to calculate 

•dis,soc ^ z (t) ~~ ^min , 
15 «n» ~ At j c 1 

•chg,soc ^ z (t) ~ z uiax 
min " T)At/C ' 

where z(t) is the pack SOC. 

In one embodiment of the present invention, the power predictive method can take into 
20 account more information than simply the cell SOC. For example, a Kalman filter can 
be used as a method to estimate all the cell SOCs in a pack. Besides giving the SOC, 
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Kalman filtering yields estimates of the uncertainty of the SOC estimate itself. A method 
of using Kalman filter to estimate SOC is described in commonly assigned U.S. Patent No. 
6,534,954, hereby incorporated by reference. 

Let the uncertainty have Gaussian distribution with standard deviation, as estimated by 
5 the Kalman filter, be denoted as a z . Then, the method yields a 95.5% confidence that the 
true SOC is within the estimate ±2a z and a 99.7% confidence that the true SOC is within 
the estimate ±3a z . 

This information can be incorporated into the estimate of maximum current based on 
SOC to have very high confidence that SOC design limits will not be violated. This is done 
1 0 as (assuming a 3a z confidence interval): 

•dis,soc = ( Z k(t) — 3cTg) — Z m j n 

max >* At/C 

,chg,soc _ ( z k(t) + 3a z ) — Z mdx 

Zmin >* ~ " v At/C 
2 Calculation Based on Voltage Limits 

Besides taking SOC limits into account, embodiments of the present invention correct a 
1 5 limitation in the prior art HPPC method for applying voltage limits (steps 12 and 22 of Figs. 
1A and IB). In the HPPC method, if the cell model of equation (1) is assumed, and that 
# chg and R dis are the cell's Ohmic resistances, then equation (2) and equation (3) predict 
the instantaneously available current, not the constant value of current that is available for 
the next At seconds. If cases where At is large, the result of the calculation poses a safety 
20 or battery-health issue, as the cells may become over/under charged. 

To overcome this problem, an embodiment of the present invention uses the following 
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cell model: 

v k (t + At) = OCV(z k (t + At)) -Rx i k (t), (7) 

This modifies the previous cell model in equation (1). Note that this model cannot be directly 
solved in closed form for the maximum current i k since z k (t + At) is itself a function of 
5 current (cf. (4)) and OCV(-) is a nonlinear relationship. Note that other cell models can be 
used as well. 

Two method embodiments are directed to solving (7) for the maximum absolute value 

of ijb(i). 



2. 1 Method I: Taylor- series expansion 

1 0 The first method uses a Taylor-series expansion to linearize the equation, so that an 
approximate value of i can be solved. It is assumed that OCV(-) is differentiate at point 
z k (t) 9 which gives the result 

OCV(z k (t + At)) = ocv 

= ocv^w)-!^— j dz 

15 +R l (z k (t),i k {ri i At)/C) 



z=z k {t) 



where the first-order residual Ri(z k (t),i k (r)iAt)/C)/\\z k (t)\\ -> 0 as i k (r}iAt)/C -> 0 
in R Note that the change in SOC over At seconds is generally small, so the following 
approximation may be applied 



OCV(z k (t + At)) m 0CV(z fc (i))-i fc ^ a0CV(2) 



dz 

20 -Rxi k 



z=z k (t) 
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which gives 



•dis,volt _ 



•chg,volt _ 



OCV(^(i))-^min ^ 



At dOCV(z) 
C dz 



+ R dis 



OCV(z k (t))-v n 



c dz 



Zk(t) 



(8) 
(9) 



In one embodiment, both the function OCV(z) and its derivative dOCV (z)/dz might be 
5 computed from some known mathematical relationship for OCV(z), (e.g., Nernst's equa- 
tion) using either analytic or numeric methods, or by a table lookup of empirical data. This 
quantity is positive for most battery electrochemistries over the entire SOC range, so the 
values computed by (8) and (9) are smaller in magnitude than those from (2) and (3) for the 
same values of R dl8 and i? chg . 

1 0 The HPPC procedure compensates for its inaccuracy by using modified values of R dl3 
and J?^ 6 , determined experimentally, that approximate the denominator terms in (8) and (9). 
This can not be accurate over the entire SOC range, however, as dOC V (z) / dz is not constant, 
particularly near extreme values of z. 

Discharge and charge currents with all limits enforced are computed as (steps 16 and 26 

15 of Figs. lAandlB) 



♦dis 
''max 



= min(w,mm<K,™fiy) (10) 
and power may be calculated using the sum of all cell powers. These are equal to the product 



15 



of the maximum allowed current and the predicted future voltage. 




n, 




(12) 



►dis 



"fE<5k«*(*+a«). 



max 




(13) 



Maximum and minimum cell and pack power limits may also be imposed in this calculation. 
Note that in all equations, OCV(s), C, <w, v m \ n9 z m ax» ^min, imax, imin, # chg , and R di3 
1 0 may be functions of temperature and other factors pertaining to the present battery pack 
operating conditions. 

2.2 Method II: Using a Comprehensive Cell Model 

The method for solving (7) presented in the previous section requires less computational 
intensity. A second method embodiment of the present invention may be used when more 
1 5 computational power is available. This second method assumes a more precise mathematical 
model of cell dynamics, which might be in a discrete-time state-space form such as the 
coupled pair of equations 




(14) 




(15) 
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where m is the discrete time sample index, the vector function of time x k [m] is called the 
"state" of the system, Uk[m] is the input to the system, which includes cell current i k [m] as a 
component, and might also include temperature, resistance, capacity and so forth, and /(•) 
and g(-) are functions chosen to model the cell dynamics. Alternate model forms, including 
5 continuous-time state-space forms, differential or difference equations might also be used. 
It is assumed that there is a method to compute SOC given the model that is implemented. 

For convenience of presentation, it is assumed that the cell model is in a discrete-time 
state- space form. Also assume that At seconds may be represented in discrete time as T 
sample intervals. Then, this model can be used to predict cell voltage At seconds into the 
1 0 future by 

v k [m + T]= g(x k [m + T],u k [m + T]), 

where x k [m + T] may be found by simulating (14) for T time samples. It is assumed that the 
input remains constant from time index m to m + T, so if temperature change (for example) 
over this interval is significant, it must be included as part of the dynamics modeled by (14) 
1 5 and not as a part of the measured input u k [m] , 

The method then uses a bisection search algorithm to find i^^ k and i^^ lt by looking 
for the i k (as a member of the u k vector) that causes equality in 

g(x k [m + T],u k [m + T]), or 
0 = g(x k [m + T},u k [m + T}) -v mia (16) 

20 to find im^jfc*. and by looking for the i k that causes equality in 

g(x k [m + T],u k [m + Tf), or 
0 = g(x k [m + T],u k [m + T]) - *w (17) 
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to find i c m ^ k • ^ special case is when the state equation (14) is linear — that is, when 



where A and B are constant matrices. The model presented in Section 3, entitled "An 
Example Cell Model," is an example where this is the case. Then, for input u k constant 
5 time m to m + T, leading to 



Most of these terms may be pre-computed without knowledge of u k in order to speed 
calculation using the bisection algorithm. 

Once the SOC-based current limits and <SS3? C m computed using (5) and (6), 
10 and the voltage-based current limits i^£ k and i c ^ k lt are computed using (16) and (17), 
overall current limits may be computed using (10) and (11) (steps 16 and 26 of Figs. 1A 
and IB). Power is then computed as 



x k [m + 1] = Ax fc [m] + Bu k [m] 




jfiZ = n p f,i^v k (t + At) 



k=l 
n 9 



= n P Jl iC min9k (x k (m + T) , u*) , 



1 5 with u k containing i m f n as its value for current, and 



P2S* = n^i^it + At) 



n 3 



with Uk containing as its value for current. 
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2.2.1 Bisection search 

To solve (16) and (17), a method to solve for a root of a nonlinear equation is required. In 
one embodiment the bisection search algorithm is used for this requirement. The bisection 
search algorithm looks for a root of f(x) (Le, a value of x such that f(x) = 0) where it is 
5 known a priori that the root lies between values x\ < root <x 2 . One way of knowing that 
a root lies in this interval is that the sign of f(xi) is different from the sign of /(x 2 ). 

Each iteration of the bisection algorithm evaluates the function at the midpoint x mi d = 
(xi + x 2 )/2. Based on the sign of the evaluation, either x\ or x 2 is replaced by x m id to 
retain different signs on }{xi) and /(x 2 ). It is evident that the uncertainty in the location 
10 of the root is halved by this algorithmic step. The bisection algorithm repeats this iteration 
until the interval between xi and x 2 , and hence the resolution of the root of f(x) is as 
small as desired. If e is the desired root resolution, then the algorithm will require at most 
|"log 2 (|x 2 - xi \/e)] iterations. The bisection method is listed in Listing 1. 

2.2.2 Finding maximum/minimum current 

1 5 To determine maximum discharge and charge current for any particular cell, bisection 
is performed on (16) and (17). Bisection is incorporated in the overall algorithm as follows. 
First, three simulations are performed to determine cell voltages At seconds into the future 
for cell current = 0, %k = i m \n, and = z max . If cell voltages are predicted to be between 
Vm\n and iWx for the maximum dis/charge rates, then these maximum rates may be used. 

20 If the cell voltages, even during rest, are outside of bounds, then set the maximum rates 
to zero. Otherwise, the true maximum rate may be found by bisecting between rate equal 
to zero and its maximum value. Bisection is performed between current limits (i m \ n , 0) or 

(0,Zmax)- 
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Algorithm 1 The bisection algorithm. 

begin {bisection algorithm} 
set xi = first search limit 
set x 2 = second search limit 

set e = desired resolution on current in bisection output 
set JMAX = maximum number of bisection iterations 

> [log 2 (^)l 
let func(-) be function in which to find root 
set the search interval dx = x 2 - xi 
if ( ( func(xi)>0) { // constraint: func(xi)<0 
dx = —dx 

X\ = X 2 

} II root is now between (x u x\ + dx), and func(a;i)<0 
fforj = 1 to JMAX { 
dx = 0.5 x dx 

Zmid = X\ + dx 

if ( func(x mi d) < 0) { xi = x raid } 

if (|<te|< e) { return(zi + 0.5 x dx) } 
} II have bisected too many times, so return best guess 
return(a;i + 0.5 x dx) 
end{bisection algorithm} 



3 An Example Cell Model 

An example cell model for the present invention power estimation methods is presented 
herein, with illustrations given to show the performance of the two methods compared to 
the prior art PNGV HPPC method. The cell model is a discrete-time state-space model 
5 of the form of (14) and (15) that applies to battery cells. The model, named "Enhanced 
Self-Correcting Cell Model," is further described in the article "Advances in EKF LiPB SOC 
Estimation," by the inventor, published in CD-ROM and presented in Proc. 20th Electric 
Vehicle Symposium (EVS20) in Long Beach CA, (November 2003) and is hereby fully 
incorporated by reference. It is understood this model is an example model only and that a 
1 0 variety of suitable alternate models can be used. 
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The "Enhanced Self-Correcting Cell Model" includes effects due to open-circuit- voltage, 
internal resistance, voltage time constants, and hysteresis. For the purpose of example, the 
parameter values are fitted to this model structure to model the dynamics of high-power 
Lithium-Ion Polymer Battery (LiPB) cells, although the structure and methods presented 
5 here are general. 

State-of-charge is captured by one state of the model. This equation is 



where AT represents the inter-sample period (in seconds), and C represents the cell capacity 
(in ampere-seconds). 

1 0 The time-constants of the cell voltage response are captured by several filter states. If 
there is let to be nj time constants, then 



The matrix Af e R n f xn f may be a diagonal matrix with real- valued entries. If so, the 
system is stable if all entries have magnitude less than one. The vector 5/6 R n f xl may 
15 simply be set to n/ "l"s. The value of n/ and the entries in the Af matrix are chosen as 
part of the system identification procedure to best fit the model parameters to measured cell 
data. 

The hysteresis level is captured by a single state 



z k [m + 1] = z k [m] - {7iiAT/C)i k [ml 



f k [m + 1] = A f f k [m] + B f i k [m]. 
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where 7 is the hysteresis rate constant, again found by system identification. 
The overall model state is 



Xk[m] = 



AH' h k [m] z k [m] 



where the prime symbol (') is the matrix/vector transpose operator. The state equation for 
5 the model is formed by combining all of the individual equations, above. Note, that at each 
time step, the state equation is linear in the input 



u k [m] = 



i k [m] 1 



which speeds the prediction operation. 

The output equation that combines the state values to predict cell voltage is 

1 0 v k [m] = OC V (z k [m]) + Gf k [m] - Ri k [m] + Mh k [m] , 

where G G R lxn f is a vector of constants that blend the time-constant states together in the 
output, R is the cell resistance (different values may be used for dis/charge), and M is the 
maximum hysteresis level. 

The open-circuit- voltage as a function of state-of-charge for example Lithium Ion Poly- 

1 5 mer Battery (LiPB) cells is plotted in Fig. 3. This is an empirical relationship found by 
cell testing. First, the cell was fully charged (constant current to 4.2V, constant voltage to 
200mA). Then, the cell was discharged at the C/25 rate until fully discharged (3.0V). The 
cell was then charged at the C/25 rate until the voltage was 4.2V. The low rates were used to 
minimize the dynamics excited in the cells. The cell voltage as a function of state of charge 

20 under discharge and under charge were averaged to compute the OCV. This has the effect of 
eliminating to the greatest extent possible the presence of hysteresis and ohmic resistance 
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in the final function. For the purpose of computations involving OCV, the final curve was 
digitized at 200 points and stored in a table. Linear interpolation is used to look up values 
in the table. 

The partial derivative of OCV with respect to SOC for these example cells is plotted in 
5 Fig. 4. This relationship was computed by first taking finite differences between points in the 
OCV plot in Fig. 3 and dividing by the distance between points (i.e., Euler's approximation 
to a derivative). The resulting data is too noisy to be of practical use, as shown in the gray 
line of Fig. 4. It was filtered using a zero-phase low-pass filter, resulting in the black line 
of Fig. 4, which may be used in the power calculation. This relationship is also digitized 

1 0 at 200 points, and linear interpolation into the table of values is used when computations 
requiring this function are performed. 

Other parameters are fit to the cell model. In particular, the model employs four low- 
pass filter states (n/ = 4), a nominal capacity of 7.5Ah, and an inter-sample interval of 
AT = Is. There is very close agreement between the cell model voltage prediction and the 

1 5 cell true voltage. This is illustrated in Fig. 5, which is a plot showing the voltage prediction 
using the cell model of the present invention. For this figure, the cell test was a sequence of 
sixteen UDDS cycles, performed at room temperature, separated by discharge pulses and 
five-minute rests, and spread over the 90% to 10% SOC range. The difference between true 
cell terminal voltage and estimated cell terminal voltage is very small (a root-mean-squared 

20 (RMS) voltage estimation error of less than 5mV). To better illustrate the model's fidelity, 
refer to the zoom on one UDDS cycle in the 50% SOC region, shown in Fig. 6. The SOC as 
a function of time is plotted in Fig. 7, which is a SOC trace for cell test. The graph shows 
that SOC increases by about 5% during each UDDS cycle, but is brought down about 10% 
during each discharge between cycles. The entire operating range for these cells (10% SOC 

25 to 90% SOC, delineated on the figure as the region between the thin dashed lines) is excited 
during the cell test. 
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4 Comparing Maximum Power Calculations 

The PNGV HPPC power estimation method gives a result that is a function of only SOC. 
Therefore, it is possible to graph available power versus SOC to summarize the algorithm 
calculations. The first method proposed (Method I: Taylor Series Expansion Method) in this 
5 patent disclosure is also possible to display in this way. Estimated power is only a function 
of SOC, dOCV/dz (also a function of SOC), and static limits on maximum current and 
power. The second method (Method II: the Comprehensive Cell Model Method), however, 
dynamically depends on all states of the system. Two systems at the same state of charge, but 
with different voltage time-constant state values or hysteresis state levels will have different 

1 0 amounts of power available. To compare power computed by the three methods, dynamic 
tests must be conducted. 

For the following results, a pack of LiPB cells is assumed to be with n s = 40 and n p = 1. 
The data to fit the models was collected from prototype hand-made cells jointly developed 
by LG Chem (Daejeon, Korea) and Compact Power Inc. (Monument, Colorado). Limits 

1 5 for the power calculations are listed in Table 1. Each cell has a nominal capacity of 7.5Ah, 
and At was ten seconds for both charge and discharge. 

Table 1: Parameters for power calculation examples. 
Parameter Minimum Maximum 



v k (t) 


3.0V 


4.35V 


h(t) 


-200A 


200A 


z k (t) 


0.1 


0.9 


Pk(t) 


— oo 


oo 



First, the PNGV HPPC method and Method I of the present invention are compared 
in Fig. 8, which is a plot comparing static maximum power calculations as functions of 
SOC for the PNGV HPPC method and Method I of the present invention. The black curves 
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correspond to charge power, and the gray curves correspond to discharge power. Note that the 
absolute value of power is plotted to avoid confusion due to sign conventions. Considering 
first the calculations of charge power, it is evident that the PNGV HPPC method produces 
similar values to Method I in the mid-SOC range. The slight differences are due to the fact 
5 that the 10-second # chg value used for the PNGV method and the derivative-modified i?* 6 
for Method I are not identical. Outside the mid-SOC range, the graph shows that Method I 
ramps power down in the neighborhood of z max to avoid over-charging the cell, whereas the 
PNGV method has no such limit. At very low SOCs, the PNGV method over-predicts how 
much power is available since there are no current limits applied to the calculation. The 
1 0 Method I estimate is automatically lower due to the large derivative in the denominator of 
the calculation. This causes an anomaly near zero SOC where the method under-predicts 
the available charge power. However, since the cell will not be operated in this range, this 
is not a concern. 

Considering now the discharge power curves, the comparison shows that Method I 
1 5 imposes limits on discharge power to ensure that the cell is not under-charged, whereas 
the PNGV method does not. In the SOC range from about 15% to 35%, the two methods 
predict similar powers. For SOC above about 35%, the power predicted by Method I 
saturates because the maximum discharge current limit of 200A has been reached. The 
PNGV method does not consider this limit. At SOC around 99% the graph again shows 
20 an anomaly in the Method I calculation where power is under-estimated due to the large 
derivative term. This apparent glitch is not a problem since the cell will not be operated in 
this range. 

Figs. 9 through 13 show how the two main voltage-limit based methods of power 
estimation of the present invention (Method I and Method II) compare to the prior art 
25 PNGV method in the dynamic cell tests shown in Fig. 5. Fig. 9 is a plot showing that 
discharge power capability estimates for cell cycle test comprising sixteen UDDS cycles 
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over an SOC range of 90% down to 10%. Fig. 10 is zoomed-in plot of Fig. 9, showing 
about one UDDS cycle. Fig. 1 1 is a plot showing charging power capability estimates for 
cell cycle test comprising sixteen UDDS cycles over an SOC range of 90% down to 10%. 
Fig. 12 is zoomed-in plot of Fig. 11, showing about one UDDS cycle. Again, the absolute 
5 value of power is plotted. 

In the discussion that follows, the results of Method II are considered to be the "true" 
capability of the cell. This assumption is justified by the fidelity of the cell model's voltage 
estimates, as supported by the data in Fig. 6. Fig. 9 shows that the three methods produce 
similar estimates. In particular, Methods I and II appear to be nearly identical when viewed 

1 0 at this scale. At high SOCs, the PNGV HPPC method predicts higher power than is actually 
available (by as much as 9.8%), and at low SOCs the PNGV HPPC method under-predicts 
the available power. Only the methods of the present invention include SOC bounds, which 
explain why their predictions are so different from the PNGV HPPC estimates at low SOC. 
If the vehicle controller were to discharge at the rates predicted by the PNGV HPPC method, 

15 the cell would be over-discharged in some cases (lowering its lifetime), and under-utilized 
in other cases. Figure 10 zooms in on Fig. 9 (same region shown as in Fig. 6) to show 
greater detail. In this region, the three methods produce nearly identical predictions. A 
notable feature of Method II, however, is that it takes into account the entire dynamics of 
the cell when making a prediction. Therefore, the strong discharges at around time 237 and 

20 267 minutes draw the cell voltage down, and allows less discharge power than the other two 
methods which only consider SOC when making their estimate. 

The three methods are also compared with respect to charge power, shown in Fig. 11. 
At this scale, the estimates appear nearly identical. Again, the PNGV HPPC method does 
not consider SOC limits, so over-predicts charge power at high SOCs. It also over-predicts 

25 power at low SOCs as it ignores the increase to charge resistance at low SOC. A zoom of 
this plot is shown in Fig. 12, which accentuates the differences between the predictions. 
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Here, it can be seen that the strong discharges at around time 237 and 267 minutes allow 
for greater charging power, as the voltage will not quickly change. 

5 Conclusion 

While the methods described herein, and the apparatus for carrying these methods into 
5 effect, constitute preferred embodiments of the present invention, it should be recognized 
that changes may be made therein without departing from the spirit or scope of the present 
invention, which is defined in the appended claims. For example, the steps 10, 12, 14 
disclosed in Fig. 1 A can be executed in different orders or used in different combinations 
and steps 20, 22, 24 disclosed in Fig. IB can be executed in different orders or used in 
1 0 different combinations. Also, various cell models can be substituted for the purpose of 
estimating the maximum absolute charge/discharge power of a battery/battery cell. 

A method and apparatus for the calculation of power capability of battery packs using 
advanced cell model predictive techniques has been described in conjunction with one or 
more specific embodiments. The invention is defined by the following claims and their full 
1 5 scope of equivalents. 
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